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Abstract 

We propose a theory for the underdoped hole-doped cuprates, focusing on the "nodal-anti-nodal 
dichotomy" observed in recent experiments. Our theory begins with an ordered antiferromagnetic 
Fermi liquid with electron and hole pockets. We argue that it is useful to consider a quantum 
transition at which the loss of antiferromagnetic order leads to a hypothetical metallic "algebraic 
charge liquid" (ACL) with pockets of charge — e and +e fermions, and an emergent U(l) gauge 
field; the instabilities of the ACL lead to the low temperature phases of the underdoped cuprates. 
The pairing instability leads to a superconductor with the strongest pairing within the — e Fermi 
pockets, a d-w&ve pairing signature for electrons, and very weak nodal-point pairing of the +e 
fermions near the Brillouin zone diagonals. The influence of an applied magnetic field is discussed 
using a proposed phase diagram as a function of field strength and doping. We describe the 
influence of gauge field and pairing fluctuations on the quantum Shubnikov-de Haas oscillations in 
the normal states induced by the field. For the finite temperature pseudogap region, our theory has 
some similarities to the phenomenological two-fluid model of — 2e bosons and +e fermions proposed 
by Geshkenbein, Ioffe, and Larkin [Phys. Rev. B 55, 3173 (1997)], which describes anomalous 
aspects of transverse transport in a magnetic field. 
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FIG. 1: Square lattice Brillouin zone map showing the electron and hole pockets of conventional 
spin density wave theor y 10 ' 11 ' 12 ' 13 of antiferromagnetic order at wavevector (jr, it). The hole pockets 
are centered at the wavevectors K„ = (±7r/2, ±7r/2) (where v = 1 ... 4), and the electron pockets 
are centered at Q a = (vr, 0), (0, ir) (where a = 1,2). The present paper describes the influence of 
quantum and thermal fluctuations in the orientation of the antiferromagnetic order on the electron 
and hole pockets. We will find strong pairing of the electron pockets at Q a , which induces a weak 
'proximity effect' pairing of the hole pockets at K„, with all the pairings consistent with a d x 2_ y 2 
pairing signature. 

I. INTRODUCTION 

A remarkable consensus has emerged in recent experiments^^ 1 ^ 1 ^ 1 ^ on the enigmatic 
underdoped region of the hole-doped cuprate superconductors. These experiments reveal 
a clear "dichotomy" between the low-lying electronic excitations near the nodal points of 
the (i-wave superconductor (i.e. near the wavevectors K„ in Fig. [1]) and the higher energy 
excitations near the "anti-nodal" points (i.e. near the wavevectors Q a in Fig. [I]) The nodal 
quasiparticles have a pairing energy which decreases with decreasing doping, and they form 
coherent quasiparticles which display characteristic interference patterns in scanning tun- 
neling microscopy (STM) observations. In contrast, the anti-nodal excitations have a larger 
gap which increases with decreasing doping, and they appear to be excitations of a state 
with STM modulations characteristic of a valence bond solid. 5 

Theoretically, a number of numerical studies^^£ of the Hubbard model have also pre- 
sented evidence for the nodal- ant i- nodal dichotomy at intermediate energy scales. These 
results are connected to ad hoc theoretical models^^'i*^ - involving "Fermi arc" and/or 
electron/hole pockets which violate the traditional Luttinger theorem on the area enclosed 
by the Fermi surfaces. A central point behind the analysis of the present paper is that 
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theories with such ad hoc violations of the Luttinger Fermi area law are fundamentally 
incomplete. Using arguments building upon the non-perturbative proof of the Luttinger 
theorem~ it was argued^ 1 ^ that metallic states with non-Luttinger Fermi surfaces must 
have "topological order," by which we mean there must be additional collective excita- 
tions associated with an emergent gauge field. Such collective excitations are crucial in the 
description of such exotic conducting states. Specific theories^ 1 ^ 1 ^ of conducting states 
with non-Luttinger areas, labeled "algebraic charge liquids" (ACL), have been provided: in 
models appropriate for the cuprates, these states were obtained across quantum transitions 
involving the loss of antiferromagnetic Neel order. Furthermore, the formalism developed 
to describe the quantum ACL is also useful for describing the "liquid" state obtained when 
the antiferromagnetic order is lost by thermal fluctuations. 

We begin presentation of our results by recalling spin- density- wave (SDW) studies of the 
onset of antiferromagnetic order in the doped cuprates.— ^ These works are expressed in 
terms of a vector SDW order parameter Ng (with I = x,y,z), measuring the spin-density- 
wave at wavevector (tt, 7r), which can perturb the Fermi surface of a weak-coupling band 
structure; we will restrict our attention here to the commensurate Neel SDW, and recent 
work by Harrison^ has shown that a similar Fermi surface structure is obtained for incom- 
mensurate SDW order. The theory for the transition from the SDW ordered state ((Ng) 7^ 0) 
to the non-magnetic state ((Ng) = 0) is expressed in terms of an effective action for space- 
time fluctuations Ng. The state with (Ng) 7^ has "small" Fermi pockets: hole pockets 
centered at the and electron pockets centered at the Q a . This conventional, ordered 
antiferromagnetic state will also be present in our theory below. In the spin-density-wave 
theory, the non-magnetic state with (Ng) = has a "large" Fermi surface which obeys the 
conventional Luttinger theory, and the transition from the small Fermi pockets state to a 
"large" Fermi surface co-incides with the loss of SDW order. 

In our theory below, the physical properties of the SDW ordered state are qualitatively 
identical to those in the spin-density-wave theory. However, we express our theory for the 
loss of SDW order not in terms of the vector Ng order parameter, but in terms of a bosonic 
spinor z a which is related to Ng by 

Ng = z>U*p, (1-1) 

where the a are the Pauli matrices. Then, the state with (z a ) 7^ is the same as the spin- 
density wave state with (Ng) 7^ 0. However, an important advantage of the formulation in 
terms of the z a is that we can describe the electron spin in terms of its components quantized 
along the direction of the local Neel order, simply by performing a SU(2) rotation defined by 
the spinor z a . This facilitates a description of the non-mag netic state^iSfi with ^ = Q> 
which is a topologically ordered ACL that retains key aspects of the "small" Fermi surface 
structure, as summarized in Fig. [2j and will be discussed in detail below. The z a formalism 
also efficiently describes the non-magnetic state obtained when SU(2) invariance is restored 
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FIG. 2: Non-superconducting ground states of our theory as a function of the coupling s which 
tunes the strength of the SDW fluctuations (see C z in Eq. (|Aip ): s is controlled by varying the 
doping 5 and by choosing different cuprate series. We expect that at zero applied magnetic field, 
H , these phases are pre-empted by superconductivity; the phase diagram as a function of H and 
5 is shown in Fig. [3J 



by thermal fluctuations. 

The primary motivation for the present paper comes from the recent experimental evi- 
dence for electron pockets in the hole doped cuprates at large magnetic fields: Shubnikov-de 
Hass (SdH) and de Hass-van Alphen (dHvA) oscillations^ 7 - 1 ^ 1 ^ 1 ^ 1 ^ indicate carriers in small 
pockets, and Hall conductivity measurements^ have been used to argue that these are elec- 
tron pockets. One of our main claims is that an algebraic charge liquid consisting of a pocket 
of charge — e fermions at the Q a wavevectors, and of charge +e fermions at the Kj, wavevec- 
tors, provides the underlying quantum state for the description of the underdoped cuprates, 
and also for the thermal fluctuations of the more classical liquid state in the 'pseudogap' 
regime. Speculations along these lines were also made in Ref. |26|. This ACL preserves the 



full symmetry of the Hamiltonian. Instabilities of this ACL involving the onset of SDW 
order, superconductivity, and charge order will be key for the description of the underdoped 
regime — see the proposed phase diagram in Fig. [3j In particular, the onset of supercon- 
ductivity removes low energy fermionic excitations which suppress monopole-instantons in 
the gauge field, and so is likely to lead to a confinement transition; this confinement physics 
has been studied in model systems in earlier work—, and we will not discuss it further here. 

We shall pay particular attention here to the pairing of the — e pocket. We note that 
a phenomenological model of pairs of electrons near the Q a wavevectors, with charge — 2e, 
was considered by Geshkenbein, Ioffe, and Larkin, 17 and we will discuss the connection to 
their model further in Section [Til As we will describe there, the charge carriers in this pocket 
experience a strong attractive interaction, which causes them to form an s-wave paired state. 
However, the resulting superconducting state actually has <i-wave pairing for the physical 
electrons,— as will be reviewed in Section [UJ Furthermore, we will show that the pairing of 
the fermions near the K„ wavevectors is very weak, and has nodal points along the Brillouin 
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FIG. 3: Proposed zero temperature phase diagram for the under doped cuprates as a function of 
doping, S, and a magnetic field H perpendicular to the layers; the insulating phases at very small 
8 are not shown. This phase diagram combines our present results with those of Refs. [10 and 
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The value of <5 C! and the scale of 5 will be differ for the various cuprate series. The superconducting 
phases are shaded and labeled dSC. The dashed line indicates the normal state phase boundary 
of Fig. [2] which is pre-empted by the onset of superconductivity. The physics of the confinement 
and possible valence bond solid (VBS) order in the dSC state is discussed in Ref. [26|. At higher 
S, there is a transformation to the physics of the "large" Fermi surface state, which we do not 
describe here. Only the ACL phase above has unconventional 'topological' order; all other phases 
have conventional order and associated excitations. While the ACL can be a stable ground state 
at large H, it is possible it is pre-empted by a conventional Fermi liquid. 



zone diagonals. Thus our theory contains a low density of strongly paired charge carriers, 
a d-wave pairing signature with nodal points, and a nodal-anti-nodal dichotomy: these are 
all attractive phenomenological features. 

We shall describe the loss of superconductivity in this state by an applied magnetic field, 
H, and the appearance of SdH oscillations at H > H c2 , using the phase diagram in Fig. 
Both the SDW metal and the ACL metal exhibit SdH oscillations. Specific predictions for 
the H and temperatures (T) dependence of these SdH oscillations will be provided, which 
can potentially be compared with experiments. Because our charge ±e fermions do not carry 
spin, a key feature of the SdH oscillations will be the absence of a linear Zeeman splitting 
in the Fermi surface areas, both in the ACL and in the phase with SDW order. 33 

Another effect of the applied H field is that it can induce or enhance the SDW order 
for the fields at which superconductivity is present — see Fig. [3j This has been discussed 
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theoretically in some detail,— and in Section [TT] we will show that our present pairing theory 
provides a specific mechanism for the competition between antiferromagnetism and super- 
conductivity This field-induced enhancement of SDW order is also clearly observed in the 
La series of hole-doped cuprates.— The samples with 5 = 0.1 display^ the be- 
havior expected for 5 < 5 C ; the samples with 5 = 0.144 have observed^^ 1 ^ the physics 
predicted^ for 5 > 5 C , including the dSC to dSC+SDW transition at a non-zero H. Thus 
for the La compounds we expect that 5 C is between the quoted values. Very recently, field 
induced SDW order characteristic of 5 > 5 C has also been reported 41 - for YBa2Cus06.45. 
The quantum oscillation experiments^ 1 ^ 1 ^ 1 ^ 1 ^ 1 - 1 ^ are on the same (or closely related) com- 
pound, and so it is plausible that they are in the metallic SDW phase of Fig. [3j as others 
have also suggested.— Our predictions for quantum oscillations and pairing instabilities will 
also extend to this metallic SDW phase. 

Our theory also offers a natural starting point for a description of the finite T "pseudogap" 
region of the underdoped cuprates. The strongly paired — e pockets lead to an effective 
description in terms of charge — 2e bosons which can exist above the superconducting critical 
temperatures.— After a duality mapping, this leads to a theory of a vortex liquid which can 
capture both the "phase" fluctuations and the possible instabilities to varieties of charge 
order.— ^ We will explore some of the thermoelectric transport properties of such a model 
in Section [TVl 

The physics of the lightly-doped Mott insulator as described by a t-J model is usually 
implicated in the description of cuprate superconductivity. 45 From this perspective, our use, 
following a proposal in Ref. 0, of charge — e carriers for the hole-doped cuprates may seem 
unacceptable. It is often argued that such excitations are only present across a Hubbard 
gap of energy U, and the limit U — > oo has been taken by a Gutzwiller projection of such 
carriers.— In response to this potential objection, we have the following responses: 

(i) We draw the reader's attention to the electron-doped cuprates, which presumably have a 
similar value of U, and for which both electron-like and hole-like carriers have been observed 
in photoemission experiments on the state with SDW order.— We are assuming here 
that similar physics applies to the hole-doped cuprates. 

(ii) While the upper Hubbard band is indeed separated by an energy of order U, a "Kondo 
resonance" of these states can be present at the Fermi level, and this is described by our 
charge — e carriers. Indeed states from both the upper and lower Hubbard bands are 
required^ 1 ^ for the eventual appearance of the Fermi-liquid Luttinger ( "large" ) Fermi sur- 
face state in the overdoped regime, and so we claim it is not surprising that both bands have 
precursors in the underdoped regime. 

(Hi) We note recent arguments by Comanac et al.,— based upon optical conductivity data 
and dynamic mean field theory arguments, that the effective U in the hole-doped cuprates 
is not as large as is commonly assumed. 

(iv) As we will review below, the electron pockets reside in a region of the Brillouin zone 
where the pairing force is strongest. Thus for these momenta, it will pay to acquire states 
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from the upper Hubbard band to benefit from the increased pairing energy. The effective 
density of carriers in the charge — e pockets will be larger than would have been assumed 
without accounting for the pairing. Indeed, even if the chemical potential is below the bot- 
tom of the band of — e carriers, the pairing interaction will induce a non-zero density of — e 
carriers. 

The outline of the remainder of our paper is as follows. In Section [TTJ, we describe our 
model for the underdoped cuprates. The primary actors are the electron pockets near Q a 
which pair strongly. We show how this model leads naturally to a <i-wave superconducting 
pairing for physical electrons, with only weak pairing and nodal excitations along the zone 
diagonals near the K„ points. Section II III will consider quantum oscillations in the normal 
state obtained by applying a strong magnetic field. We will describe the corrections to 
the Lifshitz-Kosevich formula from gauge field and pairing fluctuations. Section IIVI will 
discuss additional experimental consequences of our theory: a two-fluid model for transverse 
thermoelectric transport in the pseudogap region. 

II. THE MODEL 

The starting point of our analysis is an expressiorJ 4 -^^ for the electron annihilation 
operator ^/ Q (r) (where a =|, \ is a spin index) in terms of continuum fermionic fields F va 
and G aa , which reside in the vicinity of the wavevectors K„ and Q a respectively: 

4 2 

* a (r) = v^^Te^Fi + ^Z- g J2e lQar G aa . (2.1) 

v=l a=l 

Here Zf t9 are non-singular quasiparticle renormalization factors which depend upon mi- 
croscopic details. As will be described explicitly below, the fermions F va and G aa are in 
turn expressed in terms of a bosonic spinon field z a , and spinless fermions which carry 
the electromagnetic charge. The phases with (z a ) ^ have been shown^ to be conven- 
tional SDW-ordered states:— 1 ^ the excitation spectrum co-incides with that obtained in 
spin-wave/Hartree-Fock theory. However, the utility of the parameterizations below is that 
the same formalism can be easily extended across the quantum transition at which Neel 
order is lost, and we reach a phase with (z a ) = 0. There is substantial recent numerical 
evidence^^ 1 ^ 1 ^ that the ^-based theory correctly captures the low energy fluctuations 
across the Neel-disordered transition in insulating model systems. One of our main assump- 
tions will be that this description in terms of the z a is a valid starting point for describing 
the loss of Neel order in the doped cuprates. The theory for this Neel disordering transition 
also involves an emergent U(l) gauge field = (A T ,A), which is connected to the gauge 
field of the CP 1 model; the z a carry unit charge under A^. Here [i is a spacetime index 
extending over the spatial co-ordinates x,y and the imaginary time co-ordinate r. 

For the electronic excitations near the K. v , we need the electromagnetic charge +e "holon" 
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annihilation operators f qp , where q = ± and p = 1, 2. Here q is the "charge" under and 
p is a "valley" index; note that although there are 4 pockets near the K. v in Fig. [TJ a proper 
counting of degrees of freedom requires only two valleys. The complete expressions for the 
F va at all for K. v points in terms of the f± p are 24 

= { f h' 2 ) , ( F jf )=*.( f k' 2 ) , (2-2) 




J -1,2 J \ r 3,4| / \ J-l 



where 



n z = r T . (2.3) 

The physical content of this parameterization is simple: the ± indices of the f± p are the 
spin components quantized along the local SDW order, and these are rotated by the SU(2) 
matrix TZ Z to a fixed quantization direction by the z a . Note that pockets separated by 
the SDW ordering wavevector of (jr, it) are parameterized by the same degrees of freedom, 
and they differ only sign of the /_„ operator. Eq. (12. 2p is the same as the parametrization 
proposed in the semiclassical theory of lightly doped antiferromagnets by Shraiman and 
Siggia. 56 As discussed in previous work— in the non-SDW phase with (z a ) = 0, the 
parameterization in Eq. (12.21) and the coupling in Eq. (1A3I) lead to electron spectral functions 
which are not centered at K„; once Neel order has been disrupted, there is no special 
reason for the electronic spectrum to be pinned at the magnetic Brillouin zone boundary. 
The computed^ 4 - electron spectral functions have a "Fermi arc" structure, similar to those 
observed experimentally. An additional mechanism for Fermi arc behavior is from the phase 
fluctuations of the superconducting order, and these effects will appear in our theory from the 
"Josephson" term introduced in Eq. (12.111) between the f± p fermions and the pairs formed 
out the states near Q a ; a recent work— has examined classical thermal phase fluctuations 
present at high temperatuers, and our formulation allows for a systematic consideration of 
quantum phase fluctuations at low temperatures. 

Indeed, our primary focus here is on the electronic excitations near the Q a wavevectors. 
For these we need electromagnetic charge — e "doublon" annihilation operators g q , where 
q — ±. The g± will be the central actors in our analysis. Note that the g q do not carry any 
valley index, and the two charges of g± specify all the fermionic degrees of freedom at all 
the Q a in Fig. [H The g± are related to the physical electrons by Eq. (12. ip and 26 





(x2T 



G 



21 




K z \ y , (2.4) 



where the SU(2) rotation 1Z Z was defined in Eq. (12.31) . In the SDW state, the ± indices of 
the g± fermions (and also of the f± p fermions) become equivalent to the f, j spin indices 
quantized along the direction of the Neel order i.e. the g± are conventional electron opera- 
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tors.— However, in the phase with spin rotation invariance preserved, ± gauge charges can 
be interpreted as sublattice indices which determine the sublattice on which the fermion is 
predominantly (but not exclusively) located. 

We will carry out our analysis in the framework of an effective field theory for the g± 
coupled to the emergent gauge field. The complete Lagrangian for our field theory has 
the following structure (field theories for bosonic spinons and spinless fermions were also 
considered in early work^^^^) : 



C = C a + C z + Cf + C fg + C 



c 



+ 9! 



[d T — iA T + iea T ) — 



(d T + iA T + iea T ) — 



2m* 



2m* 



f~-zg + C'A 



fV - iA) 



^9+9-9-9+ ■ 



9+ 

9- 



(2.5) 



We have only written out explicitly the Lagrangian C g which involves the g± fermions, and 
which will be the basis for almost all the computations in the body of this paper. The term 
Cf g coupling the f± p and g± fermions will be described below; all other terms have been 
discussed previously^>2£ and are recalled in Appendix [A] In Eq. (12.51) . a T is the external 
electrostatic potential whose coupling shows that both g± carry charge — e. The fluctuations 
of a T are controlled by the action 



Sa = ^ ldT 



d 2 k 

^|k||a T (k,r)| 2 



(2.6) 



which leads to the Coulombic repulsion e 2 /r between all the g± particles. The magnetic 
dipole interactions associated with fluctuations of the electromagnetic vector potential a can 
be safely ignored. The g± carriers have any effective mass m*, and experience a chemical 
potential \x. 



A. Fermion pairings 

1. g± pairing 

We will be especially interested in the pairing of the g± fermions as described by C g . 
Indeed, we will present arguments below in favor of the proposition that an s-wave pairing 
of the g± is the primary pairing instability of the underdoped cuprates: the pairing of the 
f± p fermions, and of the physical electrons ty a will be shown to follow from it. 

Eq. (12. 5p already includes an attractive contact BCS interaction, A, between the g±. This 
attraction is permitted by the underlying symmetries,— and so can be written down on phe- 
nomenological grounds. More physically, the longitudinal component of the A^ gauge force 
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provides an important component of the attractive force between the g+ and g_ fermions: 
this is simply the attractive "Coulomb" force between two opposite charges. This will be 
Thomas- Fermi screened by the compressible fermion state to an attractive force with a range 
of order the Fermi wavelength. It is clear that this force prefers an s-wave pairing between 
the g± fermions. An additional contribution to the s-wave attractive force comes from the 
term C zg in Eq. (1A5j) . Integrating out the z a spinons, we find a contact attractive interaction 
~ -A 2 . 

However, the key source of the s-wave pairing of the g± is the force associated with the 
transverse components of the gauge field. As long as we are in the phase without SDW 
order, these remain long-range and unscreened. The nature of the fluctuations are similar 
to those of the fermionic U(l) spin liquid 62 or the Halperin-Lee-Read state.— We have the 
following propagator for the transverse part of the gauge field: 



(A(q,^)^i(-q, -w,)) 



°H q 2 



' (2-7) 



Xq 2 + l\uj\/q + A AF 



In our case, the effective gauge-field propagator contains contributions both from the z a 
spinons and from the charge carrying fermions f± p , g+. The spinon contributions to the 



susceptibility \ from Eq. ( 1A1I) were discussed in Ref. [64j. The fermions yield x = [6vr 2 z/] , 
with the effective density of states, v = Jl/tt, determined by the reduced mass of the holons 
and doublons, Ji = m*mf/(2m* + mf) (here my is related to the masses in Eq. (1A2I) . and the 
factors of 2 arise from the valley degeneracy) . The damping term 7 comes from the Landau 
damping of fermions and is given by the sum of two Fermi-momenta 7 = (v^f* + 2p^^ / (2ir). 
Finally, the "mass" term Aaf arises from the Higgs mechanism in the state with SDW order 
with 



^af~K*«)| 2 , (2-8) 



as discussed in Appendix [A] 

Pairing due to transverse gauge forces has been considered previously in the context 
of spin liquids. Because the magnetic force between two oppositely directed currents is 
repulsive, a pairing between fermions could occur only in unusual channels,—* 6 ^ in particular 
in the "Amperian" channel where the fermions on the same side of the Fermi surface pair 
up.— However, in our case note that the g + and g_ carry opposite A^ gauge charges, and 
so the magnetic force is attractive in the traditional s-wave BCS channel of pairing between 
fermions on opposite sides of the Fermi surface. Indeed, this problem of pairing by transverse 
gauge forces between Fermi surfaces of opposite charges has been considered previously 
by Bonesteel, McDonald, and Nayak, 67 and by Ussishkin and Stern 6 ^ in the context of 
double layer quantum Hall systems each at filling fraction v = 1/2. In this quantum Hall 
problem, the electrons in the two layers have opposite gauge charges with respect to an 
"antisymmetric" U(l) gauge field whose flux measures out-of-phase density fluctuations in 
the two layers, and their s-wave BCS instability leads to a paired quantum Hall state. An 
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Eliashberg analysis of such a pairing instability due to transverse gauge forces was carried 
out in these works,— and their results can be related to our problem. An important result 
obtained in these studies was that while the low-energy gauge fluctuations lead to very 
singular electron self-energies in the normal state (including non- Fermi liquid behavior), 
they are not^ 1 ^ pair-breaking; the pairing instability remains very strong. This should be 
contrasted with the behavior near ferromagnetic quantum critical points, where there is 
a similar anomalous self-energy in the normal state, but the ferromagnetic fluctuations are 
pair-breaking to p-wave superconductivity.—^^ For our problem, the estimate of the s-wave 
pairing temperature is T p0 ~ 7 2 / (m* 3 x 2 )- Using the values of x and 7 quoted below Eq. (12.71) . 
and ignoring the spinon contribution to x> we arrive at the simple estimate T p o ~ E-p, where 
Ep = p F /(2m*) is the Fermi energy for electrons. To the extent we can work within the 
context of C g in Eq. ( 12. 5p . we can understand this estimate on dimensional grounds. Note 
that in the non-magnetic phase, the only dimensional parameters appearing in Eq. (12.71) are 
associated with the Fermi surface, and there is no arbitrary coupling constant in the coupling 
between g± and A^. In this respect, this problem is similar to the three-dimensional Fermi 
gas at a Feshbach resonance. Consequently, the mean- field pairing temperature T p0 can only 
be of order the available energy scale, which are the Fermi energies. In reality, the actual 
value of T p o will be also influenced by the spinon contribution to x, the Coulomb repulsion 
e 2 /r between the g±, the contribution of the f± p to the A^ polarization, and the value of A. 

Given the quenching of the transverse gauge propagator in Eq. ( 12. 7p in the phase with 
(z a ) 7^ 0, we can expect that the pairing instability will become weaker in the SDW ordered 
state. This then sets up a natural and appealing mechanism for the suppression of T p0 after 
the onset of SDW order. Indeed, it offers a basis for the theory of "competing orders"— 
which has many attractive phenomenological features. 

We have now established that C g has a strong pairing instability to a state where 

( 9+ (k)g4-k)) = A g , (2.9) 

where we can take the pairing amplitude A g to be independent of k near the Fermi level. 
Then, what is the pairing amplitude for the physical electron operators in Eq. (12.11) ? We 
assume, for simplicity, that we are in a non-magnetic state where (z^zp) ~ 8 a p. Then from 
Eq. (12T41) we obtain 

(G la (k)G 1/3 (-k)) = -(G 2a (k)G 2 p(-k)) ~e a pA g] 

(Gi«(k)G 2/J (-k)) = 0. (2.10) 

Comparing with Fig. [H we see that this is precisely the pairing signature expected for d-wave 
pairing of the electrons; see also Fig. HI 
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FIG. 4: Sign of the d-wave pairing amplitude superimposed on the electron and hole pockets of 
the SDW metal. Note that these signs correspond to s-wave pairing of the electron pockets and 
p-wave pairing of the hole pockets, as explained in the text. 



2. f± p pairing 

Finally, we turn to the pairing of the f± p fermions. Just as was the case for the g± 
fermions, the gauge forces will prefer a s wave pairing of oppositely charged f± p fermions. 
However, there is a competing "proximity" effect arising from the paired g± fermions. This 
proximity effect is due to a Josephson coupling between g± and f± p pairs. The form of such 
a Josephson term is tightly constrained by the projective symmetry group (PSG), which was 



presented in detail in Refs. I24|j26l and the transformations needed here are listed in Table [H 
An analysis based on Table fl] shows that the simplest allowed coupling between the g± and 
f± p pairs which is invariant under the PSG is 



£fg ~ ~ iJ f9 



9+9- 



f+t D x f-i - /+i D y f-l + f +3 D x f-2 + f+2 Dy f-2 + H.c,(2.11) 



where Jf g is the Josephson coupling, Di = di — qAi is the co-variant derivative acting on a 
field with charge q, and a Di b = aDjb — (Dia)b. Note that Table fl] does not permit any 
term without a spatial gradient. From the structure of the Josephson coupling in Eq. ( 12. lip 
we see that the proximity effect induces a p-wave pairing of the f± p fermions. Thus there is 
frustration in the f± p pairing, with the gauge forces preferring s-wave. 

At the microscopic level, a computation of the pairing of the f± p fermions in the SDW 
ordered state has been carried out by Sushkov and collaborators.— iZ&Z&ZI They showed that, 
for a suitable range of parameters, the long-range spin-wave interaction preferred a p-wave 
pairing. In our formulation this long-range attraction between the f± p is mediated by the 
Shraiman-Siggia term in Eq. (lA3p . Note that this Shraiman-Siggia term does not apply to 
the g±, and so the corresponding interaction is absent there. 

On the basis of our arguments above, and the complementary microscopic 
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TABLE I: PSG transformations deduced from the PSG of the f± q in Table II in Ref. m, the PSG 
of the g± in Table III in Ref. |26I . and the PSG of A in Eq. (12) of Ref. [73|. The transformations 



are T x : translation by one lattice spacing along the x direction; '■ ^0° rotation about a dual 
lattice site on the plaquette center (x — > y,y —* 
(x y). 



/2 

); I^ ual : reflection about the dual lattice y axis 



computations^^i^i^, we conclude that the s-wave pairing of the g± fermions is the dom- 
inant instability, and the Josephson coupling in Eq. (12. lip induces a sympathetic p-wave 
pairing of the f± p fermions. A key point is that the gauge forces will be pair-breaking 
towards this p-wave pairing,— ^22. an d consequently the f± p pairing amplitude will be quite 
weak. 

Specifically, combining Eqs. (12. lip and (12. 9p . we deduce that the proximity-effect pairing 
of the f± p fermions induced by the g± fermions has the p-wave form 

(/ +1 (k)/_ 1 (-k)) ~ (k x -k y )J fg A g ; 
(/ +a (k)/_ 2 (-k)> ~ (k x + k y )J fg A g ; 

(/ +1 (k)/_ 2 (-k)> = 0, (2.12) 

where the momentum dependencies are a consequence of the spatial gradients in Eq. ( 12. lip , 
and the pairing amplitudes between fermions with like charges are zero. Finally, from 
these results and Eq. (12.21) . we can deduce the pairing of the physical electron operators in 
F va in Eq. (O): 

(F la (k)F 3/3 (-k)) ~ e aP (k x - k y )J fg A g ; 

(F 2a (k)F 4/3 (-k)> ~ e aP (k x + k y )J fg A g (2.13) 

and all other F va pairings vanish. A glance at Figs. Q] and H] shows that these are precisely 
the pairings associated with a c?-wave pairing signature of the physical electrons. The 
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momentum dependencies in Eq. (12.131) shows that the pairing amplitude changes sign across 
the Brillouin zone diagonals. Also, the vanishing of the pairing along these zone diagonals 
shows that there will be gapless "nodal" fermionic excitations. 

We close this section by noting that the structure of the Josephson coupling in Eq. 02. 11 j) 
is closely connected to that appearing in the model of Geshkenbein, Ioffe and Larkin. 1 - They 
considered a phenomenological model of charge — 2e bosons, representing pairs of electrons 
near the Q a , coupled to fermions near the K„. Their boson-fermion Josephson coupling had 
a matrix element which changed sign along the Brillouin zone diagonals. Identifying their 
boson b as b ~ g+g~, we see that Eq. (12.111) also shares these features. 

III. QUANTUM OSCILLATIONS IN THE NORMAL STATE 

This section will consider the low temperature transport properties of the paired state of 
the g± described by C g which is driven normal by a strong applied magnetic field, H. 

We begin by a simple estimate of the depairing field, if p2 (0), associated to the pair- 
ing temperature, T p o. In the absence of a complete Eliashberg theory of the influence of 
the transverse gauge fluctuations, we will be satisfied here with an estimate based on the 
weak-coupling BCS theory result for the upper critical field of a clean two-dimensional su- 
perconductor (note that the value of the upper critical field depends on the purity-^2 and 
dimensionality of the system^ 1 ^). For the purpose of numerical estimates, we ignore here 
the quantum oscillation phenomena in the transition point itself, which as discussed be- 
low may lead to a reentrant behavior. Within these assumptions, we have the following 
BCS formula^ 1 ^ that we associate with the "quantum depairing field" (here, we restore the 
fundamental physical constants): 



where e is the electron charge, c is the speed of light, m* is the effective mass of carriers in 
the pocket, ks is the Boltzmann constant, Tp = p F / \2k-Qm*) is the Fermi temperature for the 
electrons in the pocket, p F is their Fermi momentum, and 7_e = 1.781 ... is the exponential 
of Euler's constant. The onset of quantum oscillations in the cuprate experiments^*^ 1 ^ 1 ^ 1 ^ 
is identified here with the quantum depairing field, H p2 (0). The former is about H p2 (0) ~ 
50 Tesla. The quantum oscillation measurements also provide information for the effective 
electron mass, which appears to be of order free electron mass or a few times larger (the exact 
values vary in experiment) and for the area of the electron Fermi surface, which is estimated 
to be a few percent of the total area of the Brillouin zone, which in turn is determined by the 
lattice constants for YBCO. This information allows us to extract all necessary parameters. 



eff p2 (0) 




(3.1) 



m*c 



1e ft T F ' 
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The Fermi temperature is related to the frequency of quantum oscillations, FsdH, as follows: 



T F = jJ* F sdH , (3.2) 

where $o — (2.07 x 10~ 15 ) T • m 2 is the flux quantum. The first factor in Eq. (13. 2p contains 
only fundamental constants and is equal to k ™^$ = 1.33 K/T (K and T correspond to the 



units of Kelvin and Tesla respectively). Using Eq. (13.11) . we can write the following relation 
between the pairing temperature and the quantum pairing field and the Fermi temperature: 



W^felW), (3.3) 



7T 2 /Cb V m 

where /ib = eh/(2\e\c) is the Bohr's magneton. Converting all quantities to the units of 
Tesla and Kelvin relevant to the interpretation of experimental data and using the actual 
values of the corresponding physical constants, we can express the Fermi temperature for 
the electrons in the pocket and the corresponding pairing temperature as follows 

T F K « ( [F SdH T] (3.4) 



m 



and 



/n 24m 

T p0 K « Xl —^ 1 Pp K][# p2 (0) T], (3.5) 



where [T F K] is the Fermi temperature expressed in Kelvin, and [if p2 (0) T] and [-F SdH T] 
are the quantum critical field and the period of the Shubnikov-de Haas oscillations ex- 
pressed in Tesla. We emphasize that the above estimates assume the applicability of the 
weak-coupling BCS theory and a circular Fermi surface for the electron pocket. There- 
fore, Eqs. (13.11 - 13.51) are not expected to determine exactly the numerical coefficients, 
but should provide correct order-of-magnitude estimates. Using the available experimen- 



tal data, e.g. from Ref. [27J, which estimates the period of oscillations to be -FsdH ~ 530 T 
and onset of oscillations (the critical pairing field in our theory) as H p2 (0) ~ 50 T, we get 
the following relation between the zero-field pairing temperature and the Fermi tempera- 
ture: T p0 [K] « y/m e /m* a/ 12 [K] [T F K] and T F « (m e /m*) 700 K. Finally, the estimates 
for the actual numerical values of T F and T p o depend on the effective mass for electrons. 
Various experiments report different values for the latter, m* ~ (1 — 3)m e . As explained 
below in Sees. IIII Bl and IIII C\ one should be careful in extracting the effective mass from 
the temperature dependence of the oscillations in this phase, because there may be other 
effects due to superconducting and gauge fluctuations, which will change the temperature 
dependence of the amplitude in the Lifshitz-Kosevich formula. In addition, if the Fermi 
surface and/or the quasiparticle weight factor are anisotropic, it would also modify the 
effective temperature dependence in the Lifshitz-Kosevich formula^ 2 - In particular, if the 
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anisotropy of the quasiparticle renormalization Z-factor, Z p , is not taken into account, the 
effective mass extracted from the quantum oscillation measurements will overestimate the 
actual effective mass by the factor of (1/Z P ) FS ,— where the angular brackets imply aver- 
aging over the Fermi surface. However, if we now assume that the effective mass for the 
electronic excitations in the pocket is of order free electron mass (as suggested by experi- 
ments), m* ~ m e , then we get the Fermi temperature of the electron pocket Tp rs 700 K, 
the pairing temperature T p0 rs 100 K, and the zero-temperature BCS superconducting gap, 



The relatively large ratio between the electron pairing temperature and the Fermi tem- 
perature (Tpo/Tp) ~ (1/7) justifies our earlier conclusion about strong Cooper pairing in 
the electron pocket. We note further that a complete description of the finite temperature 
pseudogap region likely requires the inclusion of further interactions between the Cooper 
pairs. A particularly interesting possibility appears within an effective model where the 
paired electrons in the pockets interact on the lattice. This type of model for the electron 
pairs, (g + g_), may have a superconductor-to-insulator phase transition with a Mott-type 
gap of order, J, accompanied by development of charge order associated with the density of 
bosons 43 and the monopole Berry phases^ 3 - 1 ^ 3 - in La (see Appendix IA4I) . 

In the following, we will describe the Shubnikov-de Haas oscillations in the resistivity 
at H > H P 2(T). For reference, we recall the Lifshitz-Kosevich formula for the oscillatory 
component of the resistivity p (retaining only the lowest oscillation harmonic) 



where p\\{H = 0) is the Drude resistivity in zero field, X(T) = 2tt 2 T/uj c , t is the elastic 
scattering time from impurities, u c = eH/(m*c) is the cyclotron frequency. 

We begin our analysis of transport by discussing the fate of the Ioffe-Larkin composition 
rule in our system in Section IIII Al In Section IIIIB[ we will describe the corrections to the 
Lifshitz-Kosevich result from the fluctuations of the gauge field, while the influence of 
pairing fluctuations will be discussed in Section IIII CI 

A. Ioffe-Larkin composition rule 

Before turning to the computation of the transport properties at H > H p2 (T) in the 
subsections below, we need to discuss an important, but technical, issue. In previous studies 
of spin-charge separation in the cuprates, a crucial ingredient in the computation of the 
physical conductivity was the Ioffe-Larkin composition rule.— This states the resistivities 
of the spinons and charge carriers add to yield the physical resistivity. In our present 
situation, there is a crucial difference from the models considered in these works: our g± 
fermions carry opposite charges under the internal gauge field, and the same charge under 



A/k B = (tt/te) T p0 « 200 K. 




(3.6) 
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the electromagnetic gauge field a M . In contrast, the previous theories had ho Ions carrying 
only a single charge under the analog of A^. An immediate consequence for our theory is that 
the cross-polarization operator between the two gauge fields, H-Aa, vanishes identically: the 
g + and <?_ fermions induce opposite polarizations which cancel each other. (More formally, 
this can be seen by the PSG of the A gauge field^, which changes sign under translation 
by a lattice site, and so cannot couple linearly to a at long wavelengths.) In other words, in 
the presence of an applied electromagnetic field a M , the physical current is carried equally by 
the g + and the g-. In this current carrying state, the A^ currents of the g + and g- travel in 
opposite directions, leading to no net internal gauge current; consequently, there is also no 
spinon current. The final conclusion is then very simple: the physical conductivity is just 
the sum of the conductivities of the g + and g_, and the traditional Ioffe-Larkin rule does 
not apply to our model. 

B. Gauge field fluctuations 

The influence of gauge fluctuations on magnetotransport was examined in the context 
of the v — 1/2 quantum Hall state,— ^ and here we will adapt these earlier results to 
our problem. This analysis was carried out using the quasiclassical method, in which the 
gauge field fluctuations are treated as a random static "magnetic" field which influences the 
cyclotron motion of the fermions. We will follow the same method here. 

There are two potential sources of the random field. In the quantum Hall case, the most 
important source was the local field induced by the Chern-Simons term from the impurity 
potential. This source is absent in our case, as we do not have a Chern-Simons term. Indeed, 
the PSG of the A field 73 - shows that only impurities which locally break time-reversal can 
induce a non-zero flux of A; we will assume that such impurities are absent. An important 
consequence is that the amplitude of the SdH oscillations in Eq. (13.61) remains unaffected at 
T = by the presence of the gauge field. 

The second source of the random field was the thermal fluctuations of A. For a random 
field, h — V x A, with equal time correlations given by 

(h(r)h(r>)) = U(\r-r>\), (3.7) 

Mirlin et al— showed that the SdH oscillations in Eq. ( 13.61) are suppressed by a factor 
exp(— Sh) where 

S h = 7iR 2 c ^ ^-Jl(qR c )U(q). (3.8) 
Jo 1 

Here R c = \j2E-p/m* /uj c is the cyclotron radius of fermions at the Fermi level, and f/(q) 
is the Fourier transform of U(r). In the quasiclassical limit, the equal-time gauge field 
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correlations can be evaluated from Eq. (12.71) to yield 



^ = j^rb < 3 - 9 > 

We can also deduce from Eq. (12.71) a necessary condition for the applicability of the quasiclas- 
sical approximation: the characteristic frequency u ~ {q/l){xi 2 + Aaf) at the characteristic 
wave- vector q ~ 1/R C should be smaller than T. Now, we can insert Eq. (13.91) into Eq. (13. 8p 
and obtain 2 

S h = (i? cA /AW^) K x (i^VWx) , (3-10) 

where l\ and K\ are modified Bessel functions. In the phase without SDW order, where 
Aaf = 0, we then have 

S» = ^. (3.1D 

The value of Sh decreases monotonically into the phase with SDW order i.e. the SdH 
oscillations have a larger amplitude in the SDW state. Deep in the SDW state, where 
Aaf 3> XI ' R^i we have the limiting result for Sh (which is always smaller than the value of 
S h in Eq. (J3TID) 




2m* x Aaf ^ 

The thermal suppression of the SdH oscillations in Eqs. (I3.10H3.12]) by the factor exp(— Sh) 
will combine with the factor exp(— 2tt 2 T/uj c ) already present in Eq. (13.61) . While the T 
dependencies in the two factors are the same, they are distinguished by their B dependencies. 
In particular, we have Sh ~ T / B 2 in Eq. (13. lip , and this can serve as a characteristic 
signature of gauge field fluctuations in an algebraic charge liquid. 



C. Pairing fluctuations 

This subsection will describe the corrections to Eq. (13.61) in the context of a traditional 
fluctuating superconductivity computation built on BCS theory. We will not examine the 
interesting question of how the transverse fluctuations will modify the Cooperon operator. 
However, given the absence of pair-breaking effects in the Eliashberg computation,—^ it 
is reasonable to expect that the Cooperon will remain the same near the pair-breaking 
transition. In any case, we can also appeal to the onset of SDW order, which leads to 
Aaf > 0, to quench the gauge fluctuations. 

We begin with the Cooperon operator in the quasiclassical approximation as follows 

C(e, u\ r, r') = Q e +^{r - r')(/- £ (r - r') exp 



-2ie 



a ■ dl 



(3.13) 
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where G e {p) is the fermion Green's function in the absence of a magnetic field, and the 
latter (real magnetic field), H = V x a, is accounted for only in the gauge factor. Note 
that in our model, the i-electrons carry opposite e* = ±1 charges with respect to the 
"internal" gauge field, A, but have the same (negative) electron charge, — e, with respect to 
the external electromagnetic field. Hence, we can take advantage of the old results of Helfand 
and Werthamer,— who have proven that the Cooperon operator, C(e,u), whose kernel is 
defined via Eq. (I3.13p . is a diagonal operator in the Landau basis. Its matrix elements are 
defined simply by the expression without magnetic field, C n = (n \(e,u;q—> n)\ n), but 
with the momentum q replaced with the operator of the kinetic momentum of a Cooper pair 
7T = [q — 2zea(r)]. The corresponding matrix elements are known from the Landau problem 
in the elementary single-particle quantum mechanics, e.g., (n \tt 2 \ n) = 4eH(n + 1/2) (note 
that the Cooper pair charge is — 2e and mass is 2m*). The general expression for the 
Cooperon without a magnetic field is as follows: 

C(e, w; q) = 2th, 9 [£ {£ ~ u)] (3.14) 

J (2e - oj + (l/r)sgne) 2 + v 2 F q 2 

where v is the density of states at the Fermi level and r is the scattering time. Since we 
are interested in explaining the quantum oscillations, we assume that the latter is large and 
set it to t = oo. Note that the clean case is in fact more complicated than the disordered 
limit, because the Green's functions and the Cooperon are non-local objects (i.e., there is no 
exponential decay in space). The fluctuation propagator for superconducting fluctuations is 
an operator given by 

t{u) = 

For the purpose of describing quantum oscillations, we are interested only in the quantum 
critical point, H p2 (0), which is determined by the divergence of the matrix element at the 
lowest Landau level of the operator, Cq(0) = oo. This leads to the expression near the 
quantum pairing field as follows (here and below, the index "0" corresponds to the matrix 
element at the lowest Landau level): 

Cq(u) = -- 
v 

where r = [H — H p2 (0)] /H is the proximity to the pairing transition, and the value of the 
critical field was specified in Eq. (13. ip . 

In the expressions so far, the Cooperon dependence on the magnetic field is accounted 
for only via the gauge factor ( 13 .131) . Physically this corresponds to an approximation in 
which the motion of the Cooper pair in the magnetic field is quantized (more precisely, the 
center-of-mass motion is quantized) , but the cyclotron motion of individual electrons within 



(3.15) 



7 M 



T, 



pO. 



(3.16) 
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a Cooper pair is not accounted for. This quasiclassical approximation is valid if either 
temperature is relatively large, T 3> oj c (note that near H p2 (0), uj p2 /T p0 ~ T p0 /E F , which 
is small in the conventional weak-coupling BCS theory), or if disorder is strong enough, 
uj c t < 1, In the regime, where the oscillations are observed neither of these conditions is 
satisfied and therefore one has to take into account Landau quantization of electrons within 
a Cooper pair. This problem was considered back in the sixties, e.g., by Gruenberg and 
Giinter,— and we reiterate here the main steps to derive the oscillating transition point and 
the fluctuation effects in its vicinity. The quantity of interest is the fluctuation propagator, 
which we write as: 

AH = T , * - , (3.17) 

MH^-V — ^osc 

where £o(u) is the fluctuation propagator given by Eq. (13.16p . which does not take into 
account oscillations, and C osc is the correction to the Cooperon with the quantum oscillation 
effects, given by 

C osc = / e- eHr2 ' 4 C(v)d 2 r - C , (3.18) 



with Co is the matrix element for the Cooperon at the lowest Landau level without oscilla- 
tions and C(r) = T G{ e i r ; B)Q(—e, r; B) and the electron Green's function in a magnetic 

£ 

field is given by 

g(e,r;B) = u m&g u c £ , 6 ^,,^ A !' \ n^ , (3.19) 

frf. t£ - u c (n e + 1/2) -fi + isgn (e)/(2r) 

with L ne (z) being the Laguerre polynomial at the n e -th single- electron Landau level. It is 
these electron Landau levels that may generate de Haas oscillations above and even within 
the superconducting phase. The quantity C osc has been considered previously by Mineev^ 
and also by Larkin and one of the authors,— and it has the following form (again, we retain 
only the leading oscillation term, and drop all the higher-order harmonics) 

„ , rN 4z/ [ZJ^ 3X(T) ( 3tt \ {2<kE ¥ \ 

where X(T) = 2n 2 T/uj c and tc/(uj c t) are the familiar terms, which describe the suppression 
of quantum oscillations by the temperature and disorder (Dingle factor) correspondingly. 
Note however that there is an additional factor of 3 in these suppression terms in the 
leading oscillation harmonics for the Cooperon. This additional suppression (first pointed 
out by Mineev^ 1 ^) is due to the fact that to resolve quantum oscillations coming out of 
a Cooper pair built of two electrons, one has to "resolve" their relative cyclotron motion 
without breaking the pair. 
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Using the expression (13.201) . we obtain the fluctuation propagators follows 



1 

v 



H-H p2 (T) 



H 



C n 



-i 



(3.21) 



From Eq. (13.211) . we see that the oscillatory part in the Cooperon can be interpreted as a 
correction to the upper critical field, which too may oscillate and therefore pairing may show 
a re-entrant behavior at low temperatures. Hence in the regime where quantum oscillations 
are observed, the exact value of the "upper critical field" (even in the sense of the BCS 
pairing-depairing transition) is strictly-speaking ill-defined because there are many critical 
fields as long as oscillations are not suppressed. 

Another important circumstance has been pointed out by Champel and Mineev,— who 
argued that even below the mean- field critical field, H p2 (T), where the system is paired, one 
may see (de Haas-van Alfven) oscillations in the gapless superconductivity region, which is 
determined by the condition [H p2 (0) — H] /H <C ^/to p2 / E F In (uj p2 / E F ) in three dimensions 
and [H p2 (0) — H] / H <C ^uj p2 jE F in strictly two dimensions [here, to p2 = eH p2 (0)/ \m*c)]. 
We note that \Joj p2 / 'E F is the Ginzburg parameter, which is typically negligibly small in the 
conventional BCS systems, but is expected to be larger in the cuprates. The numerical esti- 
mates ( 13.11 - 1331) in the beginning of this section, suggest a very wide fluctuation Ginzburg 
region for the str ong ly-paired small electron pocket of our model; E.g., using the experimen- 
tal data of Ref. 27]] , we get ^uo p2 / E F ~ 1/3 for the electron pocket^. We also emphasize 
that these possible quantum oscillations in the gapless superconductivity region are different 
from the effect, which may arise from the normal vortex cores well below the critical field. In 
fact, the latter effect may be significantly suppressed in the strongly paired phase, where the 
individual vortex cores are not large enough to support a many-body electron state leading 
to quantum oscillations. 



(a) 



(b) 





(c) 
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FIG. 5: The relevant diagrams that describe the fluctuation transport near a pairing transition. 
Fig (la) shows the Aslamazov-Larkin contribution to conductivity in the clean limit; Fig (lb) 
shows the current vertex of the Aslamazov-Larkin diagram calculated in Eqs. (|3.23p and (|3.24[) : 
Fig (lc) describes the combined contribution of the Maki-Thompson and density-of-states terms. 
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We now ready to calculate the fluctuation corrections to the Lifshitz-Kosevich formula in 
Eq. (13.61) . The fluctuation conductivity in the Gaussian approximation is given by the sum 
of the Aslamazov-Larkin, density-of-states, and Maki-Thompson diagrams^ (see Fig. In 
a strong magnetic field and low temperature, all these diagrams are generally of the same 
order and in a disordered system all play an important role. However, in the super-clean 
case of uo c t — > oo, the Maki-Thompson and density-of-states terms cancel each other exactly 
at least for the non-oscillating terms 8 -^ 1 ^ 2 - and only the Aslamazov-Larkin term^ survives. 
The longitudinal part of the electromagnetic response tensor has been considered by Larkin 
and one of the authors 81 and reads 



Q{uj) = 8z/e 2 T^ 7o 2 1 (fi,^)£ (fi)i:i 



(3.22) 



where £o(ft) is the Cooper pair propagator at the lowest Landau level given by Eq. (13.211) . 
L\ is the Cooper pair propagator at the first Landau level, which near if P 2(0) is not singular 
and can be treated as a constant, and 7o i(ft,u;) is a matrix element between the first and 
the lowest Landau level of the current vertex]. The current vertex operator is defined by 
^y(Q,uj) = 7(q — > 7r; ft, u), where the latter quantity is given by the three-Green's function 
block: 

7(q;fi,w) = T Y.f (|^ v ^(p)^-(p)^- £ (q-P)- ( 3 - 23 ) 



The matrix element 701 (ft, u) was calculated in Ref. |8l| and reads 



701 (ft, w) 



y/2 1 + \u\t 



i-^r-dni + in 

2 vf 



uj\ + \uj\) 



+ 7o 



(3.24) 



where rn = l/v2eH is the magnetic length for Cooper pairs. Note that in Eq. (13.241) . 
we do not write explicitly the oscillation contribution, 7osc . The oscillations coming from 
the vertices are expected to have the usual Fermi liquid form, because the two graphs 
corresponding to the two vertices when "glued together" essentially reproduce the Drude 
conductivity diagram and therefore their contribution can estimated modulo a numerical 
coefficient in the overall amplitude. Using the technique developed in Ref. 8l| . we find the 
following main result for the fluctuation conductivity, which includes the leading oscillation 
terms as well (here, we address only the zero-temperature contribution): 



da = — 

7T 



1+77 



Posc(H) 

P(0) 



In 



1 



[H-H p2 (0)] /H-C osc / 



v 



(3.25) 



with C osc given explicitly by Eq. (13.201) 77 ~ 1 is a positive numerical coefficient, and 
p osc (H)/p(0) is the ratio of the Fermi liquid oscillation term and the Drude resistivity in 
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Eq. (13.61) . The latter does not involve any additional factors in the temperature and Dingle- 
temperature dependence as opposed to the Cooperon term, which determines oscillations in 
the transition point itself.— We observe that Eq. (I3.25P describes a decrease in the amplitude 
of the SdH oscillations upon approaching the superconducting pairing instability. 

Our result in Eq. (13.251) may have an interesting physical interpretation: The Cooper pair 
propagator C n {uj) = l/\ e g (A*(u;)A(c<;)) corresponds to the density of Cooper pairs at the 
n-th Landau level with the energy u (at least at T = 0, the Matsubara frequency can be con- 
verted into a real frequency via Feynman rotation). Near H P 2(0), only the Cooper pairs at the 
lowest Landau level play a significant role and the total density of such pairs is given by the 
integral over frequency [cf., Eq. (13T221) ]: N cp ~ / duo / '(2tt) C (u) oc - In {H/[H - H p2 (0)}}. 
Hence, the Cooper pair density scales as a logarithm in the proximity to the magnetic-field- 
tuned quantum transition and each electron within a fluctuating pair produces an oscillation 
term. If temperature and Dingle suppression are small, the functional dependence of the 
oscillatory part of the fluctuation conductivity is dominated by the Cooperon, C osc . The 
corresponding oscillation term should also survive below the pairing transition in the gapless 
superconductivity region. We note here that if such a term is detected in experiment, then 
the application of the usual Lifshitz-Kosevich formula to determine the effective carrier mass 
from the temperature dependence of the amplitude of the oscillations will overestimate the 
effective mass by the factor of three. In fact, there have been conflicting experimental reports 
about the value of the effective mass for carriers in the electron pocket and the above circum- 
stance may be relevant to this discrepancy. We expect that the temperature dependence of 
the oscillation amplitude may exhibit a crossover from 3X(T)e- 37r /^ cr 7sinh [3X(T)] at rel- 
atively low fields (when most electrons are paired, H < H p2 ) to X(T)e -7r /^ cT - ) / sinh [X(T)] 
at high fields (when most electrons are unpaired, H ^> H p2 ). 



IV. TWO-FLUID MODEL AND TRANSPORT PROPERTIES OF THE PSEUDO- 
GAP PHASE 



Our main statement is that the pseudogap phase arises from the Fermi surface recon- 
struction induced by antiferromagnetic fluctuations into a holon Fermi surface in the nodal 
region and an electron pocket in the anti-nodal region. The latter electron pocket remains 
strongly paired (but uncondensed) up to very large energy scales, which may explain why 
the underlying single-electron excitations had escaped discovery until recently. The facts 
that the small electronic Fermi surface has been "visualized" in the quantum oscillation 
experiments only in a very narrow doping range and in just one class of materials are most 
likely related to the purity of the samples. In addition, the oscillations are observed in the 
vicinity of the "magic" doping fraction, p = 1/8, where it is known that stripe and other 
competing orders are enhanced, which therefore (according to the arguments of Sec. [II]) 
should suppress the energy scales for electron pairing and hence reduce the upper critical 
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field at which the oscillations are detectable. The latter argument is consistent with the 
observation of a dip in the superconducting dome in this doping region (we reiterate that 
the actual superconducting transition temperature, T c , corresponds to Cooper pair conden- 
sation, which is different from the electron pairing temperature, T p . However, the values of 
two temperature scales are expected to correlate strongly). It is therefore likely that even 
though the electron pocket has not so far been directly observed in quantum oscillation ex- 
periments and photoemission measurements in other doping regimes and other materials, it 
does exist throughout the phase diagram of the underdoped cuprates. Its elusive nature can 
be explained by the exponential suppression of the oscillations by the Dingle factor or/and 
much larger pairing energy scales away from the magic doping fraction. 

Because the paired electron pocket is argued to be central to the phase diagram of the 
underdoped cuprates, it is important to discuss whether other unusual properties of the 
pseudogap phase (most importantly, its highly unusual transport properties) are consistent 
with the proposed mean-field state. In this section we argue that indeed many such anoma- 
lous thermal and electric transport experiments of the pseudogap phase can be qualitatively 
understood within our picture. In particular, as proposed in the early work of Geshken- 
bein et al.,— a change of sign in transverse thermoelectric response can appear naturally 
within the corresponding two-fluid model of a — 2e Cooper pair Bose-liquid and the +e hole 
Fermi liquid. Here, we discuss a complementary formulation of this model by 'dualizing' 
the bosons into vortices, while retaining the gapless fermionic excitations of the hole Fermi 
surface. Note that monopole Berry phases in La (see Appendix IA4j) are likely to play an 



important role in the vortex action,—^ particularly in the structure of any charge order 
instabilities, but we will neglect this complication here. 

Among the most unusual experiments in the pseudogap phase are the Hall and Nernst 
measurements. Here we qualitatively discuss these two effects within our two-fluid model of 
paired electrons and unpaired holon excitations. For completeness, we discuss in Appendix IB1 
the Hall conductivity of a generic anisotropic Fermi liquid. There are also corresponding 



expressions for the Nernst coefficient, as reviewed in Ref. |95[ In Appendices [C] and [Dj we 
review the dual mapping procedure and formally derive the conductivity composition rule 
of a vortex-holon two-fluid model: 



a = — - — a 

n 



'v 



l + a u (4.1) 



where ay is the dimensionless vortex conductivity matrix and <7f is the electrical conductivity 
matrix of negatively charged holons. Eqs. (14.11) has the simple physical interpretation: In 
the pseudogap phase, the electrons are paired but uncondensed, which means that Cooper- 
pair conductivity, <jcp, is finite and is given by the vortex resistivity in the dual language.— 
The total conductivity is obtained simply by adding up the Cooper pair (or vortex) and 
the holon Fermi-liquid contribution, which leads to Eq. (14.11) . Using this equation, we get 
the following expression for the Hall angle (see also, Ref. [971], where this equation was first 
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derived in a different context): 



cry + CTf Tr cf v 

where Troy = (o"v) 2 + (°v) ■ We see that both the magnitude and the sign of the Hall 
response are determined by the numerator of Eq. ( 14.21) . which involves several physically 
distinct contributions that may and most likely do have opposite signs. While the sign of 
holon Hall conductivity is determined by that of the positive hole electric charge, the sign 
of the vortex contribution is expected to be opposite at least in the strong coupling regime. 
Indeed, assuming a dilute density of electronic tightly bound pairs, we expect them to behave 
like canonical bosons with charge (— 2e), and so contribute a negative Hall conductivity - 
this is the basic picture of Geshkenbein et al. 17 

The issue of the vortex Hall angle becomes more settled in the regime away from strong 
coupling. We emphasize here that by the "vortex" transport, here we are referring to the 
contribution of the electron pockets, whose normal state transport is described in a dual 
model of vortices. Within the framework of BCS theory, vortex conductivity has been 
considered previously by Dorsey^ on the basis of the time-dependent Ginzburg-Landau 
equation, which for the case of our electronic pocket reads 

h 2 

7 (d t - 2iea T ) A g = — (V + 2zea) 2 A g + aA g - 0\ A g \ 2 A g , (4.3) 

where A g (t,r) = X e g (g+g~) is the electronic Cooper-pair wave-function, a and (3 are the 
Landau parameters, and 7 = 71 + 272 is the order-parameter relaxation time. We note 
here that the time-dependent Ginzburg-Landau theory formally has a very narrow regime of 
applicability and is expected to be quantitatively valid only in the gapless superconductivity 
region. However, it should provide a useful insight about the vortex Hall contribution in 
the crossover from strong to weak coupling. According to Dorsey,— the vortex Hall angle 
is related in a complicated way to the relaxation parameter, 7, and the structure of the 
vortex core. But generally the sign of the Hall angle (relative to that in the normal phase) is 
determined by the sign of the following parameter, (—72/71). In the clean limit of unscreened 
intervortex interactions (neutral superfluid) the situation simplifies leading to 7 = — % (which 
basically "restores" the Gross-Pitaevskii-like equation for the Cooper pair fluctuations) and 
giving the sign of the Hall angle identical to that in the normal phase, which is consistent 
with the discussion above in terms of the canonical bosons; the importance of this canonical 
boson contribution to the Hall transport was pointed out by Geshkenbein et al} 1 . In the 
weak-coupling and dirty limit, the situation is different and the sign of the Hall angle is 
determined by the Ginzburg-Landau time relaxation parameter given by (here we present 
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the result of the BCS model; see, e.g., Aronov et al..) 

^_ ! 7p91n7 £ 
1 8 2 dp. ' K ' 1 

Note that the BCS weak-coupling limit is not directly applicable to the strongly-paired elec- 
tron pocket. However, Eq. (14.41) above provides a tentative indication that upon approaching 
the regime of weak Cooper pairing (e.g., in the phase with enhanced antiferromagnetism cor- 
relations), the sign of the vortex Hall response may change. According to Eq. (14.41) . it is 
determined by the value of logarithmic derivative of the BCS pairing temperature with 
respect to the chemical potential, which is proportional to the derivative of the density of 
states at the underlying Fermi surface (in our case, the Fermi surface of the electron pocket); 
depending upon the shape of the Fermi surface, this could have a negative sign. Therefore, 
in the weak coupling limit, the signs of the Hall contributions for the h/(— 2e)-vortices and 
(+e)-holons may become the same (which could, in principle, be relevant in the vicinity of 
the "magic doping level," where competing magnetic orders are enhanced leading presum- 
ably to a suppression of electron pairing, and where the change of sign in the Hall response 
has been observed). We reiterate however that the electron pocket is expected to remain 
strongly paired in the most, possibly all, of the pseudogap phase and therefore generally 
the vortex contribution is expected to retain the electron-like sign and to compete with the 
holon contribution. 

The overall sign in the experimentally observed Hall effect will therefore be determined by 
the interplay of two physically different terms, Oy and a^, and can show reversals depending 
on the system parameters. Again, in the proximity to the "magic" doping level, one expects 
decreased disorder (i.e., decreased pinning strength) and suppression of the electron pairing 
as well, which should significantly alter the vortex contribution. We note that Taillefer et 
a/.— have reported a strong correlation between the sign-reversal in the Hall effect and the 
presence of quantum oscillations, which is qualitatively consistent with the afore-mentioned 
scenario. 

In the simplest application of the two-fluid model, the Nernst effect will be determined 
by the contribution of the h/(— 2e)-vortices and the holons. However, in reality the situation 
may be more complicated due to a strong (possibly competing) effect of the superconducting 
fluctuations^ (here we imply Aslamazov-Larkin amplitude fluctuations), which are known 
to be large compared to the Fermi-liquid terms even far above the pairing transition. Fur- 
thermore, possible charge-ordering instabilities of the vortex liquid,— and the associated 
proximity to the insulating state at p — 1/8 likely also play a role.— 
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V. CONCLUSIONS 



This paper has combined insights from recent experiments^^^^^^^^^i 1 ^ and 
earlier theoretical development 3iZiMi2£i2ii to present a fairly simple model for the underdoped 
cuprates with a firm microscopic and theoretical foundation. We start with a Fermi liquid 
state with long-range SDW order, containing electron pockets near the G a wavevectors, and 
hole pockets near the K„ wavevectors (see Fig. [p. Then we express the spin polarization of 
the electronic excitations near these pockets in terms of the local polarization of the SDW 
order; this is the content of Eqs. (I2.2p and (12.41) . The advantage of this procedure is that it 
allows to easily extend key aspects of the physics of the small pockets into the phase without 
long-range SDW order. In particular, it shows that the electronic excitations experience a 
long-range gauge force associated with an emergent U(l) gauge field A. 

A key feature of our theory is that the primary pairing instability is associated with 
electron-like pockets near the G a wavevectors. We showed that these pockets experience a 
strong attractive pairing force from the transverse gauge fluctuations, and this leads nat- 
urally to an s-wave pairing instability. However, after rotating back to the physical spin 
polarization direction via Eq. (12.41) . the resulting paired state was found to have a ci-wave 
pairing signature for the physical electrons. Next we focused attention on the Josephson 
couplings in Eq. (12. lip between the electron and hole pockets: we found that it induced 
a p-wave pairing of the holons, which was strongly frustrated by the gauge forces on the 
holons. Again, after rotating back to the physical electrons using Eq. (12.21) . this very weakly 
paired holon state was found to have a d-wave pairing signature with nodal fermionic ex- 
citations. The "nodal-anti-nodal dichotomy" is a natural consequence of this theory, with 
very different pairing physics near the G a and Ji v . 

Section II I II explored the nature of the SdH oscillations in the normal state induced by a 
strong magnetic field at low temperatures. We computed the nature of the suppression of 
these oscillations by gauge field and pairing fluctuations. 

Section [TV] explored aspects of transverse transport in the finite temperature pseudogap 
phase. Here, we also discussed connections to the boson-fermion model of Geshkenbein et 

We conclude by mentioning that to determine the correctness of our proposed theory for 
the underdoped cuprates, it would be essential to visualize the single-particle excitations 
in the electron pocket in other experiments, apart from the existing quantum oscillation 
measurements. Such new experiments should involve external perturbations, which destroy 
superconducting pairing in the electron pocket, without smearing out its small Fermi surface 
or altering the underlying magnetic or topological order that yields Fermi surface recon- 
struction. Due to the former limitation, high temperature and/or strong disorder are not 
appropriate for this purpose. Possible other means could be, e.g., to study AC transport 
in a magnetic field, looking, in particular, for cyclotron resonance effects coming from the 
single-electron excitations. Since the quantum oscillations have been observed, the mate- 
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rials are sufficiently clean to exhibit the cyclotron resonance phenomena as well. Another 
promising avenue could be to experimentally investigate non-linear transport, e.g., non- 
linear TV-curves in a magnetic field. In the vicinity of the upper critical field H p2 (T), the 
critical (depairing) current is expected to be relatively small and some manifestations of 
single-electron physics would appear at lower fields as compared to linear transport. Of 
particular interest would be also to compare the behavior of non-linear TV-curves in the 
regions with the opposite signs of the Hall conductivity. In our theory, the sign reversal in 
Hall response occurs naturally due to the competing contributions from the uncondensed 
electronic Cooper pairs and holons. This competition combined with the effects of depair- 
ing may lead a non-monotonic behavior and, possibly, sign reversals in non-linear transport 
data. 
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APPENDIX A: COMPLETE LAGRANGIAN 

Only two of the terms in the Lagrangian in Eq. (|2.5|) . C g and Cf g , were displayed in 
the body of the paper. This appendix will display the remaining terms, along with a brief 
discussion of their physical consequences. 

1. Spinons 

The Lagrangian for the z a is2^ 

C z = 1(8,- iA,)z a \ 2 + s\z a \ 2 + | (\z a \ 2 ) 2 , (Al) 

and the spinon "mass" term s tunes a transition from the SDW ordered state ((z a ) ^ 0) to 
a state with spin rotation invariance preserved {(z a ) = 0). Note that in the SDW ordered 
state, the spinon condensate induces a Higgs "mass" term, | (zq,) | 2 A 2 , for the A,. 



28 



2. Holons 

These are associated with the F va fermionic excitations near the K„ wave-vectors^ 

£/ = ££/. 

g=±p=l,2 

where j extends over x, y, mix — m 2y and m 2s = miy are the masses of the elliptical 
holon pockets, and the x and y directions are rotated by 45° from the principle square axes. 
Eq. ( 1A2I) also contains the energy Ef g , which is the analog of the "semiconductor" band 
gap, between the top of the hole (valence) band and the bottom of the electron (conduction) 
band. We expect that the value of Ef g is sensitive to the strength of the SDW order, 
decreasing as the SDW order becomes larger. 



d T — iqA T — iea T + \i — Ef g 



[dj - iqA, 



lean 



2m vl 



fqpi (A2) 



3. Spin-charge couplings 



These couple the spinons z a to the holons f± p and the doublons g±. The simplest allowed 
terms are couplings between the scalar densities, \z a \ 2 and ft p f± p , g±g±- However, there 
are also "Shraimain-Siggia" terms^ which couple operators carrying charges ±2 under A^. 
Again, the most general form of all these terms can be deduced from the PSG, as has been 
described in previous works. For the holons, the spin-charge couplings take the form 



^zf\Zot\ 



ApUv + iX *f £aP {f+if-izadxzp + f ] + 2Uz a dy Zp ] + H.C.. (A3) 



The second term is the Shraiman-Siggia term; in the SDW state, this term favors incom- 
mensurate spiral spin correlations. In the non-magnetic state with (z a ) = 0, this term moves 
the electron spectral weight away from the commensurate K„ points, to be centered on a 
"Fermi has been described in previous work.— Finally, integrating out the z a also 

leads to an attractive pairing term^ for the f± p . For the g±, the spin-charge couplings are 



t~-zg -^zgl^al ^ , 9q9q 

q 

+ A, 9 | [g\ {D x g_) z a {D~z p ) - g\ (D y g^) z a (D y z p ) 

+ e aP \gl (D y9+ ) z a * (D yZ P*) - gl (D x9+ ) z a * (D x z p *) 



(A4) 



+ H.c. 



Now the Shraiman-Siggia term has 2 spatial gradients and does not induce spiral correlations. 
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4. Gauge field 

These are induced by integrating out the matter fields, and can take a different form 
depending upon the nature of the matter excitations. When Fermi surfaces for g± are 
present, the gauge field dynamics is overdamped, as discussed in the body of the paper. 
However, in all phases, terms of the form La ~ (d^A v — d u A^) 2 are always allowed, obtained 
by integrating high energy degrees of the freedom. For the crossover to the confining state, 
we also need to consider topologically non-trivial configurations of the corresponding 
to monopole tunnelling events. These have been ignored in the present paper because the 
monopoles are suppressed by the holon Fermi surfaces, but their effects have been discussed 
earlier— in some detail. The monopoles come with Berry phases, and these are crucial 
in determining the nature of translational symmetry breaking in the confining phases. 

APPENDIX B: HALL EFFECT IN AN ANISOTROPIC FERMI LIQUID 

In this appendix, we summarize the properties of transverse thermoelectric linear response 
in a two-dimensional anisotropic Fermi liquid and present a general expression for the Hall 
coefficient. The results presented in this appendix are used in Sec. IIVI where the transverse 
thermoelectric response in the pseudogap phase is discussed. 

Consider a two-dimensional Fermi liquid with the anisotropic dispersion E = E(p). We 
also introduce the following standard notation: P ) = E(p, (fi p ) — E-p, where p = (p, P ) 
is simply the momentum in polar coordinates and Ep is the Fermi energy. This indirectly 
determines the value of the particle momentum as the function of energy and angle: p = 
p(£, P ). One can also define the particle velocity as a function of the energy and the angle: 
dE/dp a = v Q (£,0 p ). We now introduce the following identity to treat the integrals over 
momentum [below, E(p) is arbitrary function] 



where g s is a degeneracy due to an internal degree of freedom (e.g., spin in a usual electron 
liquid or a sublattice index in our theory), (...) = J Q 27r (d(j)/2Tc) ... is the average over the 
directions in the Brillouin zone (which reduces to the average over the Fermi surface if £ = 0), 
and we introduced the density of states: 



Let us also introduce the following notation for the integral that often appears in deriving 



= g, 




(Bl) 




(B2) 
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the finite temperature transport properties of a Fermi liquid: 



im\ T = I — J {£) — • (B3) 

UK hT J 4Tcosh 2 { £ /(2T)} 1 ; 

— oo 

At zero temperature, it simply gives [f{e)] T _+ = /(0). The notations defined by Eqs. fIBll) . 
(1B2I) . and (IB3I) will be used below to express the general thermoelectric response coefficients 
in a compact and intuitive form, which would allow a simple physical interpretation. 

The general theory of the Hall coefficient within the Green's function formalism has been 
considered by many author s an d we refer the read er to the corresponding literature (see, 



e.g., Altshuler and Aronov [101] and Livanov [1021 ]). In the general case of an anisotropic 
Fermi liquid with an angle-dependent scattering time, 7^ (e), one can obtain the following 
expression for the Hall conductivity in the limit of a weak magnetic field (i.e., if (uj c t) <C 1): 



d 

T l P ( £ )fe K( £ > 0p)t$(£, ^ P )f(e, P ) } 



(B4) 



In the zero temperature limit, T/E-p — > 0, and assuming an angle-independent scattering 
time, t, we get the simplified equation for the Hall response (a ^ /?): 

^ = eVi? (|K*}) FS - ( B5 ) 

The corresponding expression for the longitudinal Drude conductivity is 



e 2 r 



<^» FS • (B6) 



In the isotropic limit of a circular Fermi surface, v 2 = 2E/m, and Eq. (1B5I) reproduces the 
familiar expression for the Hall conductivity a -1 = (uj c t)ct^, with cr" = (vp/2)ue 2 T being 
the longitudinal Drude conductivity of Eq. ( 1B6I) . We emphasize here that according to 
Eq. (IB5j) . the Hall coefficient does depend on the derivative of the density of states, but the 
latter effects are not necessarily dominant. E.g., the Hall conductivity is finite even if the 
density-of-states is a constant. The sign of the Hall effect can change within a single-band 
Fermi liquid picture only if the density-of-states depends on the energy stronger than v\v 2 
in the corresponding directions, which requires a very anisotropic Fermi surface. 

APPENDIX C: DUALITY TRANSFORMATION 

The mean-field Lagrangian of this two-fluid model can be written as follows: 

C[A g , f] = C GL [A g ] + C h [f] + Cj[A g , f], (CI) 
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where A g = A c fr (g+g~) is the pairing order parameter describing electron Cooper pairs and 
£gl is the corresponding Ginzburg-Landau Lagrangian 



{-iV - 2eaf 
Am 



+ (d T — 2iea T ) 



A, + a(T)|A 9 | 2 -/3|A/ 



(C2) 



where a = (a T , a) is the physical electromagnetic field (we assumed that all effects of the 
gauge filed, A, have already been incorporated into the effective parameters) and a and (3 are 
Ginzburg-Landau parameters, which determine the modulus of the order parameter, |A S | = 
A = y/a/(2P), which we assume fixed. However, the phase, 0, of the order parameter, 
A g = A o exp(z0), is allowed to fluctuate. This leads to the XY-model for the electronic 
Cooper pairs and a Kosterlitz-Thouless transition below the BCS pairing temperature, T p0 . 
The second term in Eq. (1C1I) is the free fermion Lagrangian describing the motion of holes 



£ h [/] = r 



-iV + ea) 

2m h 



+ (d T + iea 7 



(C3) 



where the term (+ea) describes coupling of the holes to the electromagnetic field. Finally, 
the last term in Eq. (ICll) . £j[A g , /], is given by Eq. (I2.1ip and describes internal tunnelling 
between the electrons and the holes. 

The structure of the tunnelling term is given by e^ff, and it is very similar to that in 
the usual BCS mean-field model of a gapless su perco nductor (see, e.g., Ref. |97|). One can 



therefore just follow the steps used in Refs. 



97 



105 



to derive the so-called U(l) formula- 



tion^ of the vortex-fermion mixture in this model, which describes the motion of vortices 
statistically coupled to gapless fermions: The statistical interaction is that they "see each 
other" as sources of a (— 7r)-flux and therefore induce electromotive force on each other when 
moving. The first purely technical step is to introduce a new operator, h r = e l ^ r ^ 2 f r , which 
simplifies the Josephson term. The next step is to perform a duality transformation with 
respect to the bosonic Cooper-pair field, e^ r . The resulting action describing a two-fluid 
vortex-holon liquid is 



C[* v ,h]= *v 
+h* 



(-zV 



.dual 



2M V 



+ (d T + < ual + ia T ) 



2m h 



+ (0 T ~ iP T ) 



h + C 



gauge j 



(C4) 
(C5) 



where My is a vortex mass, a dual is the gauge field, which describes the Cooper pair density 
fluctuations [V x a dual = ricp(r)], and the only purpose of the fields a and j3 is to mediate 
the long-range statistical interaction between the vortices and the fermions. These fields are 



"attached" to the vortices via V x (3 = in^y^y and to the fermions via V x o; 
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This statistical interaction is described by the mutual Chern-Simons term, which is the first 
term in the following gauge part of the action 

£ g augc = -- (e,„,Aa,A/3 A ) + ±- {e^ x d v af* - 5^n Q ) 2 + ^-J c ^ vX d v af & \ (C6) 

The second term in Eq. (1C6|) describes dual gauge-field fluctuations, which physically cor- 
respond to plasmons. In a charged system, the plasmons have a gap due to the long-range 
Coulomb forces and as such the fluctuations of this gauge field are expected to be less 
pronounced than in a neutral Bose-system. Finally, the last term describes coupling to a 
physical electric Cooper pair current. 

The theory of Eqs. (lC4p and (1C6I) summarizes the following essential features of the two- 
fluid vortex-holon mixture: The vortices and the Cooper pairs "see each other" as (+2tt)- 
fluxes and induce EMFs (transverse Magnus forces) on each other when moving. Likewise, 
the vortices and the gapless fermions (holons) "see each other" as (— 7r)-fluxes and induce 
EMFs as well. We will be using this picture in derivation of the semiclassical transport 
equations below. 

A more settled issue is the question of the total effective magnetic field ("dual field") 
seen by a vortex. According, to Eqs. (JS2J) and O]) . it is = V x [ a ( dual ) - a] = 

2tt (no — |^h) . We reiterate that the Josephson term (12.111) in the action violates the indi- 
vidual conservation laws for the g- and /-particles. This means in particular that the density 
of the latter, n^, may vary depending on the phase. Another related non-trivial question is 
about the vortex statistics (with respect to each other). The direct duality transformation 
gives bosonic vortices, however other statistics are in principle possible. These are very 
interesting questions, which however are beyond the scope of the present study. Below, 
we treat the vortices semiclassically to develop a phenomenological theory of transport and 
derive the corresponding transport coefficients. 



APPENDIX D: CONDUCTIVITY COMPOSITION RULE IN THE TWO-FLUID 
VORTEX-HOLON MODEL 

The derivation of the semiclassical theory of transport in the two-fluid model is essentially 



identical to that of Ref. 97| and is based on the following equations, which describe the 
electromotive forces between the vortices, fermions, and Cooper pairs in the presence of 
currents and thermal gradients: 

j v = -a v e(j h + j C p) - A V VT; ^ D1 ^ 
jh = -o- n ejv - A h VT, 

where e is the antisymmetric tensor in two dimensions, j v , jh, and j cp are the vortex, holon, 
and Cooper pair density currents respectively and a v /h and A v are the vortex/holon di- 
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mensionless conductance and thermal conductivity matrices respectively. The latter two 
matrices generally have the form 



<7± 



a = | " | and A 

-a± a\\ 




(D2) 



The quantities of interest are the total electrical and thermal conductivity matrices for the 
system. E.g., in the absence of thermal gradients, the conductivity tensor is defined by 
jcp — £E, while the actual electric field is determined by E = ^ej v . Then, Eqs. (IDlj) can be 
easily resolved and give the following expression for the total electrical conductivity matrix: 



a 



(2e) 



a 



v + 4^ 



(D3) 



Despite the rather complicated set of arguments and transformation that have led to this 
result, the physics of Eq. (1D3j) is very simple: The total electrical conductivity in the un- 
condensed liquid phase is given by the sum of Cooper pair and hole conductivities. The 
former can be related to the vortex transport properties and due to duality is simply given 
by the vortex resistivity. One can also derive a complete phenomenological expression for 
the Peltier tensor in the two-fluid model defined via E = AVT. The exact expression for 
the Peltier tensor is rather involved, but assuming that the hole contribution to the Nernst 
effect is negligible, one can use the results of Ref. [97]. 
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